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I. INTRODUCTION 

Studies of the strong interaction at low energies are 
vital to the understanding of the structure and dynamics 
of hadrons as well as the nature of confinement. At low 
energies, the quantum chromodynamics (QCD) coupling 
becomes large and standard perturbative methods can¬ 
not be used. The main theoretical approaches at low en¬ 
ergies are lattice QCD and effective field theories, includ¬ 
ing chiral perturbation theory. Precise measurements at 
these energies provide valuable inputs and can constrain 
and test these approaches. 

The r] meson is one of the eight pseudo-Goldstone 
bosons of the broken chiral symmetry and therefore stud¬ 
ies of its decays provide a unique window into low-energy 
QCD. The r] meson is a light, neutral pseudoscalar with a 
mass of (547.862±0.018) MeV/c^ [T]. All strong and elec¬ 
tromagnetic decays of the ij are forbidden to first order, 
resulting in a relatively long lifetime and a correspond¬ 
ingly narrow width of 1.31 ± 0.05 keV. This makes the rj 
meson an ideal laboratory for the study of rare processes, 
since the suppression of many of the more abundant de¬ 
cay modes makes rare decays experimentally accessible. 

We report the measurement of the branching ratios of 
the following four rj meson decay channels: 

r] —>■ 7r^7r“7, 
r] —>■ e^e“ 7 , 
rj —>■ 7r^7r“e^e~, 
rj —>■ e^e“e“''e“ 

collected in proton-deuteron collisions at the WASA-at- 
COSY experiment using the rj —>■ decay with 

TT^ —>■ 77 (= ‘^y.y) as the normalization channel. 

Using a minimum bias data sample of rj mesons and the 
reconstruction capabilities of the WASA detector, most 
notably the charged particle tracking and particle iden¬ 
tification, we are able to isolate pure samples of several 
decay modes. It is important to note that these are the 
only current results on ry decays where the ry-mesons are 
produced in hadronic interactions, therefore they feature 
complementary experimental conditions compared to the 
results of other experiments which use photoproduction 
or e“*"e“ collisions for meson production. 

II. THE EXPERIMENT 

The WASA-at-COSY experiment was operated at the 
Cooler Synchrotron (COSY) at Forschungszentum Jiilich 
from 2006 to 2014 [5]. For the data used in this analysis, 
a proton beam with T = 1.0 GeV was impinged upon a 
deuterium pellet target. The reaction pd —>■ ^He rj is used 
to produce rj mesons at energies close to the production 
threshold, where the most favorable ratio between the 
7y production cross section and background reactions is 
found. The cross section of this production reaction is 
0.40(3)/rb [2111], meaning that up to 8 events containing 


7y mesons are produced per second at the peak luminosity 
of 2 X 10^^ cm“^s“^. 

The WASA detector is a fixed-target spectrometer, 
with a forward detector arranged to measure hadronic 
ejectiles and a central detector to detect light mesons or 
their decay products. A cross-sectional view of the de¬ 
tector appears in Fig.[2 The forward detector consists of 
an arrangement of thin and thick plastic scintillators and 
drift chambers covering the full azimuthal angle. Thick 
scintillators in the forward range hodoscope (FRH) are 
designed to measure energy loss via ionization. Thin scin¬ 
tillator layers in the forward window counter (FWC) and 
forward trigger hodoscope (FTH) provide precise timing 
information. The kinetic energy and the particle type 
can be determined from the pattern of energy deposits 
in the thin and thick scintillator layers. A proportional 
chamber system (FPC) consists of 8 layers, each with 260 
aluminized Mylar straws. Layers of the forward detector 
beyond the first layer of the FRH, including the Forward 
Range Interleaving Hodoscope (FRI) detector and the 
Forward Veto Hodoscope (FVH), were not used in this 
analysis due to the kinematics of the reaction. 

The central detector is surrounded by a CsI(Na) 
electromagnetic calorimeter with 1012 elements (SEC). 
Contained within the calorimeter is a superconducting 
solenoid providing a uniform 1 T magnetic field to the re¬ 
gion directly surrounding the interaction area. Charged 
particle tracking in this region is provided by the mini 
drift chamber (MDC), which is surrounded by an 8 mm 
thick plastic scintillator barrel (PSB) that provides pre¬ 
cise timing and particle identification. The MDC consists 
of 4, 6, and 8 mm diameter straw tubes arranged in 17 
layers that are alternatingly axial or skewed by -1-3° or 
—3° relative to the beam axis in order to provide three- 
dimensional tracking. An iron return yoke, shown in red 
in Fig. surrounds the central detector and protects 
the photomultplier tubes of the SEC from the magnetic 
field. A detailed description of the WASA detector can 
be found in Ref. Hig. 

The data for this experiment were taken over two peri¬ 
ods with four weeks in the fall of 2008 and eight weeks in 
the fall of 2009. Care was taken to provide consistent con¬ 
ditions between the two run periods. The solenoid field 
setting was 0.85 Tesla. The trigger conditions were based 
on information from the forward detector only, meaning 
the trigger was unbiased with respect to a decay mode of 
the r] meson. The trigger identified ^He ions by demand¬ 
ing large energy deposits in overlapping azimuthal sectors 
of the scintillator layers. In the case of the pd —>■ ^Her; 
reaction the ^He stops in the first layer of FRH, so the 
trigger included a veto on the signals from the second 
layer. Since fusion to a ^He represents only about 1% of 
the total cross section at this energy, the above conditions 
were sufficient to bring the trigger rate down to below a 
thousand events per second, which were recorded. 3 x 10^ 
events containing ry mesons were collected in total, with 
1 X 10^ being collected during the first period in 2008. 
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FIG. 1. (Color online) A cross-sectional scheme of the WASA detector with the beam coming from the left. Hadronic ejectiles 
are measured with the forward detector on the right while meson decay products are measured with the central detector on 
the left. Components are described in the text. 


III. EVENT SELECTION 

All decay channels are analyzed using a common analy¬ 
sis chain and settings up to the point of channel selection 
and kinematic fitting. The first step is to identify the ^He 
ion. Forward-scattered tracks are reconstructed by using 
hit patterns in the FPC and matching them to signals in 
the scintillator layers. To separate ^He ions from protons, 
deuterons, and charged 7r-mesons, the energy deposited 
in the FWC is correlated with the energy deposited in the 
stopping layer. Once the ^He ion is identified, the missing 
mass, MM(^He), can be calculated by determining the 
invariant mass remaining when the measured ^He four 
vector is subtracted from the known initial conditions of 
the beam and target. The resulting distribution is shown 
in Fig. with a peak at the mass of the 77 meson. This 
peak is composed of all decay modes of the 77 , since this 
stage of the analysis does not include any condition on 
the central detector. This initial sample contains 3 x 10^ 
events with 77 mesons. 

Particles from 77 meson decays are measured in the cen¬ 
tral detector. Tracks reconstructed in the MDC are ex¬ 
trapolated to the PSB and to the calorimeter. Clusters 
in the calorimeter that are not correlated with the tracks 
are treated as electromagnetic showers caused by pho¬ 
tons. A threshold is placed on a cluster energy of 20 
MeV to filter out low-energy background signals. The 
PSB time signals with a time resolution of 1 ns are used 
to provide a start signal for the drift-time measurement 
in the MDC. All tracks are required to pass closer than 
1 cm from the beam axis. 


xio^ 



FIG. 2. (Color online) Initial spectrum of the ®He missing 
mass. The prominent peak is due to 77 meson production, 
while the broad background distribution is from production 
of two and three rr-mesons. The dashed (red) line shows a 
fit of the background (incorporating the shapes of double and 
triple TT meson production, as described in the text) and the 
inset shows the peak after background subtraction with the 77 
mass marked by the dotted (blue) line. 


Event candidates for the 77 decay channels have to in¬ 
clude a ^He ion in time coincidence with at least the min¬ 
imum number of tracks (charged particles) and neutral 
clusters (photons) for a selected decay mode. 

Energy and momentum conservation is imposed by re¬ 
quiring all event candidates to pass through a kinematic 
fitting routine for a specific 77 decay channel. The kine¬ 
matic fit takes into account reconstruction uncertainties 
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for the different particle types as a function of angles and 
energies. A condition on the mass of the decaying rj or 
7 r° meson is not imposed. Events with a fit probability 
less than 0.1 are rejected. For certain channels additional 
conditions are applied, which are described in the later 
sections. 


IV. SIMULATIONS 

Efficiencies used for acceptance correction are calcu¬ 
lated using Monte Carlo simulations. Kinematic event 
generation is performed with the PLUTO-f-l- software 
package [5]. This contains realistic physical descriptions 
of all relevant channels. The angular distribution of 
the produced ry mesons measured in Ref. [7] for the 
pd —>■ ^He?? production reaction is used. 

For r] —>■ 7 r+ 7 r“ 7 r°, the Dalitz plot parameters from 
the Crystal Barrel [ 8 ] measurement are used. The sim¬ 
ulation of channels r/ —>■ 7 r+ 7 r“ 7 *^*^ is based on calcula¬ 
tions from Ref. [5]. The decays rj are sim¬ 

ulated using form factors calculated assuming the vec¬ 
tor meson dominance model with the transition form 
factor F[q^) = 1/(1 — where q^ is the invari¬ 

ant mass squared of the electron pair and parameter 
br, = 1.78 GeV^ [10]. 

Detector simulations are performed using the WASA 
Monto Carlo package, WMC, which is based on GEANT3 
m- Temporal, spatial, and energy resolution of the de¬ 
tector elements is implemented in the WMG using data 
to fine-tune the parameters. For example, the —>■ 77 
and ry —77 decays are used to determine the energy res¬ 
olution of the calorimeter by analyzing the two photon 
invariant mass distributions. Additionally, inactive de¬ 
tector channels are continuously monitored and mapped 
in the simulations. The number of inactive channels is 
typically less than a few percent. 

The WASA-at-COSY experiment uses an internal tar¬ 
get with frozen pellets injected at rates of several thou¬ 
sand per second [12] • Though vacuum pumps are po¬ 
sitioned as closely as possible to the interaction region, 
a certain amount of residual gas is present in the re¬ 
gion around the target which comes from the evapora¬ 
tion of pellets. This is quantified in data by selecting 
pd —> ^He 7 r+ 7 r“ events and reconstructing the vertex 
from the tt+tt” tracks. The resulting spectrum has a 
large spike in the target region (with dimensions deter¬ 
mined by the profiles of the beam and the pellet stream) 
as well as tails along the beam axis due to beam-gas in¬ 
teractions. Monte Garlo simulations of these “rest gas” 
events were performed by including the shape of the ver¬ 
tex distribution in the z direction deduced from the ex¬ 
perimental data. Over 90% of all events occur within one 
centimeter of the center of the interaction region. 


V. PARTICLE IDENTIFICATION 

Three of the final states being studied contain electrons 
and positrons in the final state. In WASA-at-COSY, 
above the kinematic threshold for double 7r°7r° and tt+tt” 
production the dominant background contributions stem 
from TT meson production. In the case of rj —>■ e“'"e “7 and 
ry —>■ e“'"e“e"'"e“, the final states are mimicked by more 
abundant channels containing charged 7r-mesons. In the 
case of ?y —)■ 7r+7r“e+e“, pions have to be identified in or¬ 
der to reconstruct the kinematics of the final state. For 
the two purposes (7r-meson rejection vs. identification) 
two slightly different algorithms are used. 

The WASA detector provides an independent measure¬ 
ment of the momentum of a charged particle in the MDC, 
as well as energy loss in the plastic scintillator barrel and 
in the electromagnetic calorimeter. Fig. [^ shows the dis¬ 
tributions of the deposited energy versus the momentum 
times the particle’s charge for energies measured in both 
the PSB and the calorimeter. The energy loss in the 
plastic scintillator is corrected for the track length of the 
particle in the detector element. In order to illustrate the 
discriminating power of this method, these spectra were 
created for events with four charged particle tracks, since 
a large number of these events contain electrons. The 
bands from charged 7r-mesons, electrons, and positrons 
are labeled. 

To utilize this information, a simple Bayesian approach 
has been developed which allows the discrimination of tt- 
mesons from using all pieces of information simul¬ 
taneously m- Two algorithms exist: a rejection algo¬ 
rithm, which considers two particles at a time and aims 
at rejecting 7r-meson pairs, and a classification algorithm, 
which considers four particles at a time and aims at as¬ 
signing particle types. 

For the rejection algorithm, pairs of oppositely-charged 
particles are considered: a pair could be either an e“''e“ 
or a 7r“'"7r“ with equal a priori probability. The proba¬ 
bility that a single particle is a 7r-meson or a lepton is 
determined from the momentum and energy losses and 
added to the posteriors using Bayes’ equation. After con¬ 
sidering both particles, the configuration with the high¬ 
est probability is chosen. A distribution of the posterior 
probabilities in data and simulation is compared in Fig.]^ 
where a clear separation between particle types can be 
seen. The distribution is made for events passing the 
pd —>■ ^Hee+e “7 kinematic fit hypothesis with a (^He) 
missing mass within ±1 MeV/c^ of the actual ry meson 
mass, in order to enhance the electron contribution. 

The graphical identification bands shown in Fig. [^ are 
represented as probabilities using neural networks from 
the ROOT TMultilayerPerceptron class |I4j . The 
neural networks are trained using simulated 7r-mesons 
and electrons tracks with isotropic directions and a 
flat energy distribution as signal and an uncorrelated, 
randomly-generated data set as background. The likeli¬ 
hood function that is used for Bayes’ equation is statis¬ 
tically determined from an independent set of simulated 
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FIG. 3. (Color online) Charge-signed momentum versus en¬ 
ergy deposit in the a) plastic scintillator barrel (PSB) and 
b) calorimeter (SEC). The bands corresponding to and 

TT-mesons are labeled. 


data. The resolution and position of the identification 
bands shown in Fig. was tuned in simulations to de¬ 
scribe the data before the neural networks were trained. 

For the ly —>■ Tr+Tr^e+e” analysis a classification algo¬ 
rithm is used, as described in Section [Xl] 


VI. PHOTON CONVERSION SUPPRESSION 

Final state electron-positron pairs originate from vir¬ 
tual photons where the corresponding radiative decay has 
a branching ratio that is about two orders of magnitude 
larger. The external conversion of real photons is sup¬ 
pressed by the design of the WASA detector, which uses a 
thin beryllium beam pipe in the interaction region. How¬ 
ever, there is still a 1 % chance that a photon will convert 
in the beam pipe or the inner layers of the drift cham¬ 
ber, producing an electron-positron pair. This means the 
magnitude of the conversion background is similar to that 
of the signal. 

In order to suppress this background, the electron- 



FIG. 4. (Color online) Posterior probabilities for events pass¬ 
ing the ^ e^e “7 kinematic fit hypothesis showing the dis¬ 
crimination of TT-mesons from electrons. The red and blue 
lines are simulations of ?7 —^ e''"e “7 and rj —> 7r"'"7r“7 respec¬ 
tively and the shaded area is the sum of both simulations. 


positron vertex position is determined from the recon¬ 
structed MDC tracks. For events where the electron- 
positron pair originates from the actual 77 meson decay, 
the vertex is close to the center of the interaction region. 
For events where the particles are the result of photon 
conversion in the detector material, the vertex distance 
will be at least equal to the radius of the beam pipe. The 
vertex distance in the plane perpendicular to the COSY 
beam, py, is represented on the y-axis in Fig. [^for sim¬ 
ulations and data. For simulations the channel 77 —>■ 77 
is shown in red while the colored spectrum represents 
77 —>■ e“''e“ 7 . The 77 —)■ 77 events are clustered starting 
at pv of about 30 mm, which corresponds to the radius 
of the beam pipe. The 77 —>■ e“'"e “7 events are clustered 
around zero. 

The invariant mass of the e“'"e“ pair at the beam pipe 
location provides additional discriminating power, as in¬ 
spired by a similar condition used in Ref. m- Normally, 
the direction of the momentum vector of a particle is set 
to be tangent to the track helix at the closest approach 
to the origin. For this calculation the momentum vec¬ 
tor is recalculated at the point where the helix crosses 
the beam pipe. This does not change the magnitude of 
the momentum vector, which is determined by the ra¬ 
dius of the helix, or the polar angle, which is determined 
by the pitch of the helix, but changes the azimuthal an¬ 
gle. The expectation is that the momentum vectors of 
an electron-positron pair originating from the beam pipe 
will be parallel at this point and will cause a peak in the 
invariant mass distribution around 2me- Particles origi¬ 
nating from a decay at the origin will have an offset in the 
azimuthal angle, causing an offset in the invariant mass 
distribution. This can be seen in Fig. where the con¬ 
version events have a distribution starting at zero while 
the 77 —>■ e~^e~'y events are offset. 
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FIG. 5. (Color online) Condition on conversion suppression 
shown for a) simulations and b) data. For the simulations, 
the multi-colored plot is for the 77 —>■ e^e “7 signal while the 
red dots are for 77 —>■ 77 where a photon converts in the beam 
pipe. The demarcation lines shown in each plot are the same. 


The selection condition used is illustrated in Fig. 

As part of a consistency check for each studied decay 
channel the selection condition was varied systematically 
and the effect on the final result was determined to be 
negligible. This involved changing both the slope and y- 
intercept of the diagonal component of the demarcation 
line illustrated in Fig. as well as the height of the 
horizontal component. 

VII. NORMALIZATION: 77 

The channel 77 —?► 7 r+ 7 r“ 7 r° has a branching ratio of 
0.2292 ± 0.0028, making it the most probable decay of 
the rj meson containing charged particles in the final state 
[T] . Due to the large branching ratio, a data sample with 


high statistics and low background could be extracted. 
The decay kinematics of this channel have been the sub¬ 
ject of detailed studies using the WASA-at-COSY 2008 
data. The results have been recently reported [TB]. In 
the present analysis the decay serves as a normalization 
channel for the less abundant processes. 

The MM(^He) distribution for events passing the kine¬ 
matic fit condition for the pd —>■ ^He 7 r+ 7 r “77 hypothe¬ 
sis is shown in Fig. Due to the large signal to back¬ 
ground ratio, no additional selection conditions are re¬ 
quired. The 77 peak after subtraction of the continuous 
background is expected to contain only a 2 % contribu¬ 
tion from 77 —^ 7 r"'' 7 r “7 with one spurious neutral cluster. 
The smooth background under the peak in Fig. is com¬ 
posed mostly of pd ^He 7 r+ 7 r“ 7 r° events with a small 
contribution from pd —>■ ^He 7 r“'" 7 r“ events. In order to 
fit the background, the MM(^He) spectra of these two 
processes are determined from Monte Carlo simulations 
assuming a homogenous phase space distribution. They 
are multiplied by a fourth-order polynomial to model the 
acceptance and possible deviations from the phase space 
distributions. An additional parameter controls the rel¬ 
ative scaling of the two- and three-pion continuum spec¬ 
tra. When the fit is performed the region ±3 (t around 
the peak is excluded. The number of signal events is then 
determined by subtracting the background function from 
the experimental spectrum in the signal area ±3tT around 
the 77 peak and integrating the resulting spectrum. Using 
this method (136, 240 ± dlOgiot) and (290, 810 ± 590stot) 
77 —>• 7 r+ 7 r“ 7 r°..j, events are found for the 2008 and 2009 
data sets, respectively. 



FIG. 6 . (Color online) Missing mass of ^He for events pass¬ 
ing the kinematic fit probability condition for the hypothesis 
pd —>■ ^He 7 r^ 7 r“ 77 . The background fit is described by a 
dashed line and is derived using the method described in the 
text. 


The same method of background subtraction is used 
for the other channels, with the exception of 77 —>■ 
e“'"e“e“'"e“ where the selection conditions are stringent 
enough to reject nearly the entire continuum background. 



















The systematic error of the fitting procedure is deter¬ 
mined by varying the method used. One method is iden¬ 
tical to that described above but includes the line shape 
of the rj peak, determined from simulations, along with 
an extra scaling parameter. Third and fourth order poly¬ 
nomials are also used to model the background, with both 
the fit range and exclusion range systematically varied. 

VIII. rj —>■ e^e“7 

The ?7 —)■ e+e“7 Dalitz decay proceeds via a real and a 
virtual photon intermediate state with the virtual photon 
converting into an e“'"e“ pair. According to the vector- 
meson dominance model, the virtual photon can mix with 
neutral vector mesons. This mixing is dominated by the 
p-meson (m = 775.26 ±0.25 MeV, T = 149.1 ±0.8 MeV) 
with contributions from the tails of the w-meson (m = 
782.65 ± 0.12 MeV, T = 8.49 ± 0.08 MeV) and </)-meson 
(to = 1019.461 ± 0.019 MeV, T = 4.266 ± 0.031 MeV) 
distributions. The squared four-momentum of the virtual 
photon corresponds to the squared invariant mass of the 
e“'"e“ pair, the invariant mass distribution of e"'"e“ pairs 
is affected by this mixing and the transition form factor 
can be extracted (see for example Refs. [IIIIIH]). In 
this publication we present only results on the branching 
ratio while assuming the transition form factor according 
to the vector meson dominance model. 

The branching ratio of this channel given by Ref. [1] 
is (6.9 ± 0.4) X 10“^ based on the measurements in Refs. 
[HiiinMi]- The largest data samples to date consist of 
(1345 ± 59) and (2.2 x 10^) events [T71[T8]. This channel 
was also measured in the pd —>■ ^Hep reaction by the 
CELSIUS/WASA collaboration with (435 ± 35) events 
collected [H]. 

The first selection condition is a threshold of 100 MeV 
for the energy of the photon. The cut does not reduce sig¬ 
nal efficiency, but significantly reduces the contribution 
of the pion background where tt+tt” pairs are combined 
with a spurious neutral cluster. A neutral low energy 
cluster in the calorimeter can come from noise or coin¬ 
cidental background that has not been rejected by the 
standard analysis. 

After the cut the fraction of the events in the MM(^He) 
rj peak is 70% the signal channel, 5% p —>■ 7 r+ 7 r“ 7 , and 
24% rj —>■ 77 , with a small remainder coming from rj —>■ 
7 r“'' 7 r“ 7 r°. These numbers are determined by studies of 
Monte Carlo simulations. 

In order to reduce the photon conversion background 
from the ry —>■ 77 decay, the conversion condition intro¬ 
duced in Section |VI| is applied. Simulations show that 
this reduces the contribution of ry —>■ 77 to nearly zero, 
while reducing the number of signal events by 20%. The 
background from rj 7 r+ 7 r “7 is rejected by using the 
particle identification rejection algorithm presented in 
Section |V] This reduces the background contribution to 
about 1 % with about a 10 % decrease in the number of 
the signal events. The final MM(^He) ry peak consists of 


over 98% 7 y —)■ e+e 7 events. 



FIG. 7. (Color online) ^He missing mass distribution for 
events passing all selection conditions for 7 —>■ e''"e“ 7 . The 
background fit is described by the dashed line. 


After all conditions are applied the number of events is 
determined from the rj peak content of the MM(^He) dis¬ 
tribution (Fig.[^. After subtracting the small remaining 
background from competing ry decay channels the peak 
contains (14,040 ± 120stat) V —>■ e+e “7 events. 

The systematic error was determined by varying the 
selection conditions and checking if the specific choice 
for any condition has a systematic effect on the result. 
e“'"e“ identification: The default condition is a probabil¬ 
ity of at least 50% that the particles are e+e”. This 
is varied from 30% to 70%. 

Photon conversion: Both the slope of the diagonal line 
and the height of the horizontal line composing the 
selection demarcation (see Fig. are varied sys¬ 
tematically. 

Kinematic fit probability: The condition on the probabil¬ 
ity of the kinematic fit is varied from its default 
of 10% up to 35% for the signal and normalization 
channels simultaneously. Conditions below 10% are 
not useful because the probability distribution is 
not flat in this region and includes a large amount 
of background, making extraction of the signal dif¬ 
ficult. 

Instantaneous luminosity: The instantaneous luminosity 
is monitored by recording the rates of the elastic 
scattering trigger and the pellet target every sec¬ 
ond. The branching ratio is extracted in bins of 
the elastic scattering trigger normalized to the tar¬ 
get rate. The normalization to the target rate is 
needed to account for the duty factor of the pellet 
target. 

Continuous background subtraction: Several methods of 
background subtraction are applied as described in 
Section IVIII 
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Calibration and luminosity. The data sets from 2008 and 
2009 differ only by calibration and luminosity pro¬ 
file. These parameters are included in the detector 
simulation and their systematical uncertainty is es¬ 
timated by separate analysis of the two data sets. 
This leads to assignment of 4% systematic uncer¬ 
tainty due to deviations between the two data sets. 

The methodology from Ref. [32] was applied to com¬ 
pare the significance of a proposed systematic effect to 
the differences in statisical error of the subsets of data 
used to derive this significance. Using this method on 
the the first four conditions above, none exhibit a signifi¬ 
cance exceeding 3cr and a systematic error is not assigned. 
The final test, differences in calibration and luminosity 
profile between the 2008 and 2009 datasets, shows a 4% 
deviation common to all channels. This error has been 
included in the final results. 

It is not possible to check the effect of the background 
subtraction using this method. Therefore, the systematic 
error on the background subtraction was determined by 
performing several different fits of the background and 
taking the standard deviation. The error on the accep¬ 
tance due to the uncertainties on the transition form fac¬ 
tor is determined to be negligible by performing the anal¬ 
ysis using simulated data with several different values for 
the transition form factor. 

The branching ratio relative to 77 —>■ 7 r+ 7 r“ 7 r° is ob¬ 
tained after correcting for the respective backgrounds 
and the final acceptance of 12 . 0 % for the signal channel: 

r(?7 -)■ e+e“7)/r(77 7r+7r“7r°.^) = 

(2.97± ± 0.13,j,,) X 

The resulting branching ratio is in reasonable agree¬ 
ment with other experimental values [ 1 ] , with a precision 
limited by the systematic error resulting from the lumi¬ 
nosity profile. 

IX. r] —^ e''"e“e^e“ 

The decay rj —>■ e+e“e+e“ is closely related to the 
decay rj —)■ e+e “7 above and proceeds via two virtual 
photons. The additional electromagnetic coupling sup¬ 
presses the branching ratio of the decay by two orders 
of magnitude compared to ry —>■ e“'"e“7. The only mea¬ 
surement where this decay is observed was performed by 
the KLOE Collaboration. The branching ratio was de¬ 
termined to be (2.4±0.2±0.1) x 10“® based on (362±29) 
events [23] . 

Event candidates with at least two positively and two 
negatively charged particles measured in the WASA cen¬ 
tral detector are passed through a kinematic fitting rou¬ 
tine with the pd —)■ ^He e+e“e+e“ hypothesis. Only 
events fulfilling energy and momentum conservation at 
greater than 10 % probability are further considered. 

To suppress the charged-particle background from tt^- 
mesons the particle identification rejection algorithm in¬ 


troduced in Section [V] is used. Since the algorithm con¬ 
siders one positively and one negatively charged particle 
pair at once, an e“''e“e“''e“ event candidate is accepted if 
both pairs in the two possible combinations of oppositely 
charged particles passes the selection condition. 

Background from photon conversion in the beam pipe, 
predominantly from the reaction rj —>• e“'"e“ 7 , is sup¬ 
pressed using the method presented in Section [VT] Again 
an event is accepted if both pairs in the two possible 
combinations pass the condition. 



FIG. 8 . (Color online) The MM(®He) distribution for events 
passing all selection conditions for the rj —^ e^e“e^e“ sig¬ 
nal. The background coming from yr-meson production is de¬ 
scribed by the horizontal red line. The region between the 
vertical dashed lines was excluded from the fit of the back¬ 
ground. 

After applying these criteria, the remaining number 
of 77 events is extracted from the MM(^He) spectrum 
shown in Fig. [^ Due to the limited statistics a simpli¬ 
fied method with constant continuous background term 
is used in the fit. The range from 0.535 to 0.560 GeV/c^ 
is excluded from the fit (marked by the two dashed lines 
in the distribution). The number of remaining 77 events is 
19.7±4.9stat, which is determined by counting the events 
in the signal region after subtracting the background fit. 

Background from other 77 meson decays is determined 
from Monte Carlo simulations. The only channel found 
to contribute with at least one event is 77 —>■ e+e“ 7 . Af¬ 
ter subtraction of the 77 -decay background, 18.4 ± 4.9stot 
events remain. 

Each selection condition was studied to identify possi¬ 
ble systematic effects on the branching ratio. The checks 
for systematic effects include e+e” identification, pho¬ 
ton conversion and kinematic fit probability as described 
in Section [VlII| None of these checks produces an effect 
with a significance exceeding 3cr, so a systematic error is 
not assigned. The systematic error included in the result 
comes from the 4% error assigned due to differences in 
the calibration and luminosity profiles of the two data 
taking periods, as well as the error determined by using 
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different background fits. 

The branching ratio relative to 77 —>■ 7 r+ 7 r“ 7 r° is ob¬ 
tained after correcting for the respective backgrounds 
and the final 3.3% acceptance for the signal channel: 

r(77 —)• e''‘e“e''‘e“)/r(77 —)• 7r+7r“7r°.^) = 

{lA±0Astat±0-2sys) X 10-4 

This is only the second analysis of this channel to reach 
a finite value of the branching ratio. The result is com¬ 
patible within errors to the previous analysis |23j . 

X. T] —>■ TT^Tf-y 

The decays t] —> 7 r''' 7 r -7 and r] —> Tr+Tr-e+e- are 
driven by the same underlying mechanism, corresponding 
to anomalous terms in the QCD action. These anomalies 
are described by the Wess-Zumino-Witten Lagrangian, 
which contains two terms pertinent for the rj decays 
[SUES]. The so-called “triangle” and “box” anomalies 
describe respectively the coupling of a pseudoscalar to 
two vectors and the coupling of a pseudoscalar to two 
pseudoscalars and one vector. The names are inspired 
by the shapes of the corresponding Feynman diagrams. 
The rj Tr+Tf-y^*^ reaction is described at the lowest or¬ 
der of the chiral perturbation theory entirely by the box- 
anomaly. However, within the framework of the vector- 
meson dominance model, the triangle anomaly will dom¬ 
inantly contribute since the 7 r"'' 7 r- pair in P-wave comes 
from the meson contribution. 

Various theoretical approaches attempt to determine 
the relative contributions from these diagrams and in 
particular to predict the contribution of the box dia¬ 
gram for the two observables: the branching ratio and 
the shape of the Tr'^7r~ invariant mass spectrum liEBl- 
[32]. The channel is the second most probable rj decay 
channel to charged particles with a branching ratio of 
(4.22 ± 0.08) X 10-2 d]. It was studied by few experi¬ 
ments [53HSZ]. The two most recent results: from WASA- 
at-COSY (using the 2008 pd data) [35] and from KLOE 
m, provide the 7 r“'" 7 r invariant mass spectrum with suf¬ 
ficient precision to see an influence of the box diagram 
contribution. The branching ratios normalized to the 
rj —>■ 7 r+ 7 r- 7 r° decay from CLEO |35| and KLOE |37| col¬ 
laborations are significantly below previous values. 

Events are selected with at least two oppositely- 
charged particles and one neutral particle fulfilling the 
kinematic fit requirement for the pd —>■ ^HeTr+Tr-y 
hypothesis. At this point the content of the r] peak 
in the MM(^He) spectrum is composed of 70% rj —>■ 
TT+TT-y, with the remaining background mostly due to 
77 —>■ 7 r“'' 7 r- 7 r° events where one photon is not detected. 
This contribution can be reduced by placing a condi¬ 
tion on the missing mass squared of ^He, tt’*', and tt-, 
MM 2 (^He 7 r+ 7 r-). For the signal channel ry —)■ Tr+Tf-y 
the MM 2 (^He 7 r+ 7 r-) distribution peaks at zero, while 


for rj —)■ 7 r+ 7 r- 7 r° it peaks at the squared mass of the 7 r°. 
Rejection of the events with MM2(^He tt+tt-) > 0.005 
GeV^/c^ increases the signal content of the 77 peak to 
91%. The effect of the cut and the impact of the remain¬ 
ing 77 —>■ 7 r+ 7 r- 7 r° contribution is illustrated in Fig. 
The experimental points correspond to the 77 peak con¬ 
tent determined from the MM(^He) distributions for each 
MM 2 (^He 7 r+ 7 r-) bin. 



FIG. 9. (Color online) The experimental MM 2 (®He 7 r"'" 7 r-) 
distribution for the pd —> ^HeTy events after the kinematic fit 
probability cut for the pd —^ ^He 7 r"'" 7 r-y hypothesis (points) 
compared to MC of the ry —>■ 7 r''' 7 r-y (light shaded) signal and 
the 77 — >■ tt ’*"77-77° background (dark shaded). The red curve 
is the sum of the simulations. Kinematic variables used to 
obtain the distributions are not corrected by the kinematic 
fit. 


The non-resonant background comes predominantly 
from the pd —>■ ^He 7 r“'" 7 r- reaction where a spurious pho¬ 
ton is detected. Reduction of this contribution decreases 
both systematic error on the background ht and the sta¬ 
tistical error of the final result. A major source of the 
spurious photons comes from so-called hadronic splitoffs. 
This happens when an interaction or a decay of one of 
the charged 7r-mesons creates a secondary particle which 
leaves a signal in an isolated calorimeter module. In this 
case a spurious neutral cluster is reconstructed. 

A condition to reduce the contribution of the splitoffs is 
applied for the photon candidates with low energy which 
are close to the expected impact point of the charged pion 
track in the calorimeter. The condition was optimized to 
minimize the statistical error of the extracted number of 
the signal events. 

The MM(^He) distribution for the final selection is 
shown in Fig. |10[ The background fit is performed using 
the methods described in Section lyTTl and after subtrac¬ 
tion of the 77 —>■ 7 r’'' 7 r- 7 r° contribution to the 77 peak the 
number of signal events is (139, 760 ± 430). 

No variation exceeding Scr for the checks due to kine¬ 
matic fit probability, instantaneous luminosity, described 
in Section [VIII[ is observed. The calibration and lumi- 













11 



KLOE collaboration with 1555±52 events leading to a 
branching ratio of (2.68 ± O-OQ^jat ± O.OTgys) x 10“^ [15]. 

The channel is also interesting due to searches for a 
possible CP-violation mechanism outside of the Standard 
Model (ini HI]. It has been shown that a contribution to 
the decay amplitude from the CP-violating electric tran¬ 
sition would result in a linear polarization to the virtual 
photon. A non-zero polarization of the virtual photon 
contributes to an asymmetry of the distribution of the 
angle, (j) (the dihedral angle), between the electron and 
TT-meson decay planes in the rj meson rest frame [40] . 
The (f) angle is shown in Fig. 11 The asymmetry, is 
defined as: 


FIG. 10. (Color online) The MM(®He) distribution for events 
passing selection conditions for rj —>■ yr+Tr'y. The function 
used to subtract the continuum background is shown as a 
dashed line. 


nosity comparison of the 2008 and 2009 data sets leads 
to assignment of 4% systematic u ncert ainty as in the 
r] —)• e+e “7 analysis (see Section VIII). The continu¬ 


ous background subtraction is also investigated as for 
r] —> e"'"e“ 7 , with both polynomials and the methods 
from Section IVTll 

The two specific conditions are investigated separately: 
Missing mass squared cut: The cut was varied in steps 
in the region between the mass squared and the 
signal peak at zero. 

Splitoffs: The selection condition used to reject splitoffs 
was removed from the analysis chain and the result 
remains consistent. 

The above two tests show that the conditions do not in¬ 
troduce systematic deviations and therefore overall sys¬ 
tematic error is determined by the background subtrac¬ 
tion and the difference between the two data sets. 

The branching ratio normalized to ry —>■ 


r(77 




7)/r(?7 ■ 

0.206 ± 0.003,t,i/;,t± 0.008 sys 


TT+TT = 


XI. g —>■ 7 r"'" 7 r e''"e 

The decay g —>■ 7 r+ 7 r“e+e“ is closely related to 77 —>■ 
7 r+ 7 r “7 and corresponds to the conversion of the virtual 
photon leading to about a factor of a suppression. There¬ 
fore the measurement of this branching ratio provides ad¬ 
ditional information on the mechanism contributing to 
the parent process, g —>■ 7 r+ 7 r“ 7 *. However, the small 
decay probability ( 0 ( 10 “"^)) has made the channel diffi¬ 
cult to detect until recently. The process has been ob¬ 
served by several experiments [m [m |38l |39], but the 
only measurement of the branching ratio with statisti¬ 
cal significance more than Scr is a recent result from the 


^ V(sin (j) cos cj) > 0) — N (sin (j) cos ^ < 0) 

N{sin (j) cos (j) > 0) + N(sin (j) cos (j) < 0) 

where N{...) is the number of the decays fulfilling the 
corresponding condition. 

The theoretical upper limit for Aip is determined by 
constraints on the strong CP-violation from neutron elec¬ 
tric dipole moment measurements to be about 1% [40] . 
A previous measurement from the KLOE collaboration 
of Afj, = (— 0.6 ± 2.5stai ± l- 8 sys) X 10 “^ constrains the 
asymmetry \A^\ to be less than a few percent [15] . 



FIG. 11. (Golor online) Definition of the dihedral angle (j) for 
the g —> 7r"'"7r“e^e“ decay in the g meson rest frame. 


The analysis follows the steps outlined in Section jlllj 
using event candidates with at least two positive and two 
negative reconstructed tracks in the MDC. A unique as¬ 
pect of this decay channel is that the final state contains 
both charged 7 r-mesons and leptons. The kinematic fit¬ 
ting assumes the pd —>■ ^He 7 r+ 7 r“e“''e“ hypothesis and 
all four possible mass assignments are tested. The events 
with probability above 0.1 for at least one of the combi¬ 
nations are accepted for the further analysis. 

All four combinations for the selected events are eval¬ 
uated according to the particle identification routine de¬ 
scribed in Section]^ Additional information about decay 
angles between the oppositely-charged pairs is included in 
the algorithm. The angle between the leptons is expected 
to be small compared to the angle between the 7r-mesons. 
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This feature was previously used in Ref. |38j . The sim¬ 
ulations of the decay ry —>■ 7r+7r“e^e“ with the matrix 
element from Ref. [5] are used to determine the proba¬ 
bilities for the correct identification of the tt+tt”, 

and TT^e^ pairs as a function of the opening angle. The 
angular information is added to the probabilities, again 
using Bayes’ equation, and the conhguration with the 
highest probability is accepted. This method has been 
tested with simulations and the correct configuration is 
selected in over 90% of events. 

A significant background comes from photon conver¬ 
sion in the reactions ry —>■ 7r+7r“7 and ry —>■ 7r“'"7r“7r°. 


The conversion suppression introduced in Section VIII 


reduces the contribution of these channels to 5% of the rj 
peak in the MM(^He) distribution. The largest remain¬ 
ing background is from the ?y —> 7 r"'' 7 r“ [tt*^ —> 
decay chain, and constitutes 15% of the peak. 


e 7 ] 


using a method presented from Ref. |15j . The sin (j) cos (j) 
distribution for the selected data sample is shown in 
Fig. and compared to a Monte Carlo simulation as¬ 
suming a flat (/) distribution. The data are divided into 
sin (j) cos (j) > 0 and sin (j) cos (/) < 0 sub-samples leading to 
the two MM(^He) distributions. Due to the low magni¬ 
tude of the continuous background, the fit of the multi- 
pion background distribution uses a variety third and 
fourth order polynomials where the fit range and peak 
exclusion range are changed systematically. The rj peak 
content is obtained as the average value for the fits. The 
number of signal events in each class is obtained by in¬ 
tegrating the peak after background subtraction, further 
subtracting background from other ry decay channels de¬ 
termined from simulations, and correcting the result for 
acceptance. 



FIG. 12. (Color online) The MM(®He) distribution for events 
passing the 7 —>■ 7 r''' 7 r“e^e“ selection criteria. The back¬ 
ground fit is shown as a dashed line (red) and the shape of 
the peak from simulations of the pd —>■ ®Hery events is shown 
as a dotted line (blue). 

The missing mass distribution for events passing all 
selection conditions is shown in Fig.[^ After subtraction 
of all background channels there are (251 ± 17) signal 
events found in the combined data set. 

Due to the high statistical error all of the systematic 
effects from the kinematic fit probability and photon con¬ 
version selection conditions were determined to be neg¬ 
ligible. The 4% error due to differences in the 2008 and 
2009 data periods was nevertheless included as deter¬ 
mined from the higher statistics decays. 

The systematic error on the final result is the same 
magnitude as the statistical error. 

r(7y -J> 7r+7r“e+e“)/r(?y 7r+7r“7r°.^) = 

(1.2±0.Riat±0.Ry,) X 10-3. 

The angle between the e^e” and tt+tt- decay planes 
was determined for each event in the final event sample 



FIG. 13. (Color online) sin cj) cos (j) distribution for the data 
and Monte Carlo simulation of the rj —>■ 7 r"'" 7 r-e^e- with a 
fiat (j) distribution. 


The same sources of systematic error were considered 
as for the branching ratio analysis. Only the error on the 
continuous background fit is included in the final results 
and it is determined by the standard deviation of the 
various fits. All other errors were insignificant compared 
to the statistical error. 

The final result for the asymmetry: 

= (-1.1 ± 6.6stat ± 0.2,J,,) X 10-2. 


XII. CONCLUSION 

The obtained results on the relative branching ratios 
relative to the normalization channel ry —>■ Tr~^Tr~Tr^^ 
are summarized in Table HI The deduced value for 
F( 7 y 7 r+ 7 r- 7 )/r(?y ^ 7 r+ 7 r- 7 r°) is 0.206±0.003stat//it± 
O.OOSsys. It is in good agreement with the older exper¬ 
iments [33l [34] but is 2.6 and 2.5 standard deviations 
above the recent values from CLEO [35] and KLOE [37] 
respectively. 
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Channel Branching Ratio 

w.r.t. rj —>■ 

r] —>• 7 r+ 7 r “7 0.206 ± 0.003stat/fit ± O.OOSsys 

r, ^ e+e -7 (2.97± 0.03,t,t//it ± O.lS.ys) x 

77 — >• TT + TT'e+e" (1.2 ± o.lstat ± 0.1sj,s) X 10~^ 

77 e+e-e+e" (1.4 ± O.W ± 0.2,^,) x 10"^ 

TABLE I. Summary of experimental results for branching ra¬ 
tios relative to the normalization channel 77 —>■ 7 r^ 7 r“ 7 r°.^. 


The measured relative branching ratios can be trans¬ 
lated to absolute branching ratios by using known world 
averages from Ref. [T] for the branching ratios of 77 —>■ 
7 r“'" 7 r“ 7 r*^ and —>■ 77 . The results are presented in Ta¬ 
ble d 

The branching ratio for ry —> e^e “7 is consistent with 
the most recent Particle Data Group fit (6.9±0.4) x 10“^ 
but it is more precise by 20%. The absolute branching 
ratios for 77 —>■ 7r+7r“e+e“ and 77 —>■ e“'"e“e“'"e“ decays are 
in good agreement with the values reported by KLOE. 

[mun] 

The measured dihedral angle asymmetry, for 77 —>■ 
7 r“'" 7 r“e“'"e“ has been determined to be consistent with 
zero: = (- 1.1 ± 6 . 6 stai ± 0 . 2 ,,j,s) x 10 “^. 


Channel Branching Ratio 

77 ^ 77 + 77-7 (4.67± ± 0.19,y,) x 10'^ 

77 ^ e+e -7 (6.72 ± 0.07,t,t//it ± O.SRas) x IQ-^ 

77 ^ 77 + 77 -e+e- (2.7 ± 0.2stat ± OAsys) X 10-^ 

77 —>■ e+e-e+e-_(3.2 ± O.Qstat ± O.Ssys) x IQ-^ 

TABLE 11. Summary of experimental results for the absolute 
branching ratios, extrapolated from the relative branching ra¬ 
tio for each channel with respect to 77 —>■ 77 + 77 -[ 77 ° —>■ 77 ] us¬ 
ing the branching ratios from Ref. [I]: BR{ri —>■ 77 + 77 77°) = 
(2.292±0.028) x IQ-^ and BR( 77 ° ^ 77 ) = (98.823±0.034) x 
10 -^ 


After the collection of data presented here, WASA- 
at-COSY has collected a high statistics data sample of 
77 mesons using the proton-proton production reaction. 
This new data set is particularly important for rare de¬ 
cay studies since an order of magnitude increase in the 
number of 77 meson decay events is expected. The back¬ 
ground to signal ratio and the detector resolution are 
comparable to the presented pd data. 
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